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Raman light scattering of the two-dimensional quantum spin system SrCu2(BOs)2 shows a rich 
structure in the magnetic excitation spectrum, including several well-defined bound state modes at 
low temperature, and a scattering continuum and quasielastic light scattering contributions at high 
temperature. The key to the understanding of the unique features of SrCu2(B03)2 is the presence 
of strong interactions between well-localized triplet excitations in the network of orthogonal spin 
dimers realized in this compound. 
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Low-dimensional quantum spin systems with a quan- 
tum disordered ground state and a finite spin gap form a 
subject which is of both fundamental and applied interest 
to the physics community. This is due to the fascinat- 
ing and diverse physics of the spin-liquid state, allowing 
to address issues related to the nature of quasiparticle 
excitations and the role of strong interaction effects. It 
is known that a gapped singlet ground state is realized 
in one dimension (ID) in dimerized or frustrated spin 
chains, and in even-leg spin ladders. Of even greater in- 
terest is the existence of a gap in two-dimensional (2D) 
spin systems, because of its potential relevance for the 
description of high-temperature superconductivity. Un- 
fortunately, only very few 2D systems, such as CaV^Og, 
have been investigated, in which the spin gap is relatively 
large and not caused by anisotropies. Recently, the com- 
pound SrCu2 (B03)2 with a layered structure was identi- 
fied as a 2D S=l/2 Hcisenberg system with a unique ex- 
change topology |Q,U leading to an exact dimer ground 
state, thus providing the opportunity to study the ex- 
citation spectrum of such a model quantum many-body 
system. 

SrCu2 (1303)2 has a tetragonal unit cell (D2d) with 
Cu 2+ ions that carry a localized spin S=l/2. The spins 
form dimers which consist of neighboring pairs of planar 
rectangular CUO4 plaquettes with a Cu-Cu distance of 
2.9 A. The strength of the antiferromagnetic intradimer 
exchange coupling is estimated to be Ji = 100 K ||. The 
spin dimers are connected orthogonally by a triangular 
BO3 unit. The distance between next-nearest Cu neigh- 
bors is 5.1 A, and the interdimer exchange J2«68 K ||. 
The orthogonal arrangement of dimers thus represents a 
2D frustrated quantum spin system fjl§|. A sketch of 
the Cu-dimers is shown in the inset of Fig. 1 . The ground 
state in this exchange topology depends critically on the 
ratio J2/J1 For small J2 the system consists of 



isolated dimers and the ground state is a product of sin- 
glets, while for small Ji the model can be mapped on 
a 2D square lattice of spins and a Neel-ordered state is 
expected. It has been shown that the critical ratio of 
coupling constants that separates a gapfull and a gap- 
less state is J2/J1 = 0.7 The experimental ratio of 
J2/J1 0.68 is just below this value, placing this mate- 
rial in the dimerized phase, close to the Neel boundary. 
Thermodynamic measurements support the existence of 
a dimer ground state and a spin gap of A =34 K [QJ|] . The 
magnetic specific heat and susceptibility both show max- 
ima (T^ lax =15 K, Tmax=8 K) and a rapid decrease to- 
ward lower temperatures with exponential tails at T — > 0. 

We have investigated the magnetic excitation spectrum 
of single crystals of SrCu2 (603)2 by Raman scattering 
experiments. At low temperatures (T<cA) we observe 
well-defined low-energy singlet modes which can be in- 
terpreted as collective bound state excitations of two 
and three elementary triplets. We demonstrate that the 
appearance of the collective modes reflects the strong 
triplet-triplet interactions present in the system. The 
spectral weight evolution for intermediate (Tf»A) and 
high (T>A) temperatures can be described by the tem- 
perature dependence of the magnetic susceptibility and 
the magnetic specific heat. 

The samples were investigated in quasi-backscattering 
geometry with light polarizations in the afe-plane of the 
freshly cleaved crystal. The (a'b') axes are rotated by 45° 
within this plane. The experiments used the A = 488-nm 
excitation line of an Ar + ion laser and a laser fluence be- 
low 20 W/cm 2 . The scattered light was analyzed using 
a XY-Dilor Raman spectrometer and a back-illuminated 
CCD detector. Measurements in a magnetic field were 
performed in 90° scattering geometry. In the analysis of 
our results the main emphasis has been on the role of in- 
teractions and their influence on the magnetic excitation 
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spectrum. A detailed analysis of the phonon spectrum 
will be given elsewhere. 

In the low energy region with Raman shifts compara- 
ble to the triplet gap A=24cm _1 , drastic changes and 
a large shift of spectral weight occur with decreasing 
temperature. This evolution leads finally to the appear- 
ance of several new modes. In Fig.l Raman spectra of 
SrCu2 (1303)2 are shown to illustrate these effects at dif- 
ferent temperatures and in two scattering configurations 
with light polarizations within the a6-plane of the crys- 
tal. The spectra in the upper panel (a) present data in 
(ab):B 2 and the lower panel (b) shows data in (a'6'):Bi 
scattering configuration. Symmetry components are de- 
noted with respect to the T)2d point group using b'=a+b. 

At low temperatures (T<cA) four well-defined modes 
with energies E s — 30, 46, 56 and 70 cm' 1 appear. The 
dominant intensity of these modes is observed in the 
(a'b') scattering configuration. They neither split nor 
shift in an applied magnetic field up to 6 T as shown in 
the inset of Fig. 1(b) and therefore are assigned to spin 
singlets. The only effect of the magnetic field is observed 
as a shift of a weak intensity shoulder at 24 cm' 1 toward 
lower frequencies. This signal corresponds to the elemen- 
tary spin gap. At higher temperatures (T« A) all modes 
are strongly damped. They are replaced by a contin- 
uum of scattering with a center of gravity near 50 cm -1 
(Fig.l(a),T=7 K), corresponding roughly to the energy 
2A. For even higher temperatures (T>A) quasielastic 
scattering with a Lorentzian spectral function is detected. 
The latter two scattering intensities are observed in the 
(ah) scattering configuration. 

To understand these dramatic changes of the Raman 
spectra with temperature we compare them with observa- 
tions related to the triplet excitation spectrum and then 
map the temperature dependence of the scattering inten- 
sity onto the corresponding thermodynamic data, i.e. the 
magnetic susceptibility and the magnetic part of the spe- 
cific heat. To gain a deeper insight into the structure of 
the spectrum we will also identify the interactions leading 
to the formation of singlet bound states in the relevant 
Heisenberg model at T=0, and present estimates for their 
binding energies. 

Recent ESR jjj and neutron scattering investigations 
H on SrCu2 (1303)2 at low temperatures observed triplet 
excitations with a spin gap of A=34 K. This triplet 
branch has a very small dispersion of only 2 K point- 
ing to extremely localized excitations 0]. In addition, a 
second triplet branch A'=55 K with a larger dispersion 
of 17 K was observed. This branch can be interpreted 
as a triplet bound state of two elementary triplets (see 
discussion below). Frustration due to the interdimer cou- 
pling J2 can lead to the reduction of the ratio A'/A=1.62 
below 2 (corresponding to non-interacting magnons) j^] . 

The four modes that we observe in Raman scatter- 
ing in SrCu2 (1303)2 for T<§cA are clearly related in en- 
ergy and temperature scale to the spin gap of the com- 



pound. A phonon origin of these modes can be ruled out 
due to their temperature dependence and the nonexis- 
tence of a structural phase transition below room tem- 
perature. On the other hand, in low-dimensional spin 
systems with strong triplet-triplet interactions, well de- 
fined modes can appear below the scattering continuum 
g-jll|. In this case the light scattering process is better 
described as spin conserving (AS=0) scattering on singlet 
bound states. Such states composed of two elementary 
triplets have been observed, e.g., in the low-temperature 
dimerized phases of CuGe0 3 and NaV 2 5 PJ-pJ- 

To illustrate the mechanism of bound state formation 
at T=0 in the 2D Heisenberg model H = JyS^-Sj 
with Jij > defined in the inset of Fig. 1(a) (the Shastry- 
Sutherland model Jl5]| ), we have derived the effective 
bosonic representation in terms of triplets (tj), excited 
above the singlets, formed on the stronger (Ji) bonds: 

H = J i I>1 • *< + v X>! x • (*5 x *i) + 

« (hj) 

£ {iT (tj xuntU-tU + h - c -)} + E< f - 0>> 

i£A ieB 

(1) 

where V = —T = — Ja/2. The site indices run over the 
square lattice formed by the dimers (singlet pairs of spins 
connected by J\ bonds), and stands for nearest 

neighbors, while sub-lattice A(B) is formed by the verti- 
cal (horizontal) dimers. An important feature of Eq.(l) is 
the absence of quantum fluctuations (i.e. tj-tjj terms), re- 
flecting the fact that the ground state is an exact product 
of singlets. The triplet excitations would have been com- 
pletely localized in the absence of the T term in Eq.(l), 
and hopping appears only in order (T/Ji) e ~ {Jz/Jxf, 
leading to a small bandwidth compared to the gap HQ . 

Next we present estimates for the energies of the collec- 
tive, two-particle excitations. The zero-momentum two- 
magnon bound state with S=0, constructed by exciting 
two triplets, has the form: \^> s ) = £\ . i>f tj t\ ■ t]|0). 
Starting from the limit J2/J1 <C 1, the dominant contri- 
bution to binding comes from the two-particle interaction 
(V term in Eq.(l)), with corrections of order (J2/J1) 2 
and higher. This interaction provides an effective at- 
traction between two triplets. Assuming the triplets are 
localized the binding problem can be solved exactly, and 
the singlet binding energy defined as e s = 2 A — E% , is 
e s = Z 2 J2- Here E% is the energy of the singlet bound 
state, and we have estimated the renormalization factor 
Z by evaluating the appropriate lowest order diagrams to 
be Z as 0.75 at J2/J1 = 0.65. We have found, by compar- 
ing neutron scattering data for the one-particle disper- 
sion from with the dispersion obtained by high-order 
perturbative expansions 0] that the ratio J%jJ\ = 0.65 
is consistent with the experimental results. Indeed, for 
this ratio the theoretical bandwidth (BW=difference be- 
tween the energy at k = (%, 0) and (0,0)) and gap are: 
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BW « 0.04 Ji, A » 0.36Ji [Q, consistent with the mea- 
sured values -B14 7 « 0.3 meV and A w 3 meV ||. From 
the same analysis we estimate J2 ~ 61K w 42 cm" 1 . 
For the higher value of J2/J1 — 0.678 the bandwidth 
is approximately twice as large Q which would make 
this ratio inconsistent with the neutron scattering data. 
Putting everything together we have the estimate for the 
singlet binding energy e 5 = Z 2 J 2 ~ 24 cm" 1 . Our ex- 
perimental value for the lowest singlet mode in Fig. 1(b) 
with £f =30 cm- 1 is t s Raman « 18 cm" 1 . 

The analysis based on the localized picture some- 
what overestimates the binding since corrections of order 
(J2/J1) 2 and higher have been neglected. Such high- 
order corrections lead to the development of a strong 
dispersion as well as create several two-particle singlet 
modes with different energies, classified according to the 
group D2d- Detailed calculations have recently been per- 
formed 16 I?]] , providing quantitatively accurate results 



for the energies in the different symmetry sectors. Only 
the modes with T^{xy) symmetry contribute to the Ra- 
man intensity, and, as T — > 0, their contribution grows 
for the (a'b 1 ) and vanishes for the (ab) geometry [[l?]], con- 
sistent with our observations. In addition, the calculated 
energies jlTj are in excellent agreement with the two low- 
est (30 cm" 1 and 46 cm" 1 ) singlet modes observed in our 
experiment, which we therefore interpret as two-particle 
bound states. 

In SrCu2(B03)2 the four modes in Fig. 1(b) show ra- 
tios E s I A=1.25 - 2.9. Two of the well-pronounced 
peaks (at 56 cm" 1 and 70 cm" 1 ) are above the two- 
particle threshold (2A). The smallness of the one-particle 
bandwidth makes it possible to resolve such higher en- 
ergy peaks, since the weight of the scattering back- 
ground is expected to be small. The pronounced sharp- 
ness of the peaks indicates that they can be interpreted 
as three-particle excitations, whose existence is indeed 
possible due to the localized nature of the states. A 
three-particle singlet bound state can be constructed as: 

= Ei, jk xt j) • ^1°)' and its ener §y on thc 

same level of approximation as discussed above (localized 
triplets) is: E% = 3A — Z 3 J 2 ~ 54 cm -1 , in very good 
agreement with the energy of the mode at 56 cm -1 . The 
appearance of additional levels (as well as finite disper- 
sion) at higher order is also expected. The total Raman 
scattering intensity, reflecting the presence of two- and 
three- particle singlets with spectral weights I2 and ^3, 
is: I B {u) = h S(uj - £f) + I 3 5{u - E§). In order to 
estimate the ratio of the two intensities we notice that 
while light couples directly to the singlet two-magnon 
state, the three-magnon state has to be created via the 
action of the T (magnon number non-conserving) term in 
Eq.(l) (i.e. this term provides a vertex correction to the 
Raman operator R ~ t\ • t| ) . Consequently, we estimate 
at T=0, I3/I2 oc (J2/J1) 2 ~ 0.4, in good agreement with 
the experimental ratio of 0.32 (at T=1.5 K). Notice that 



in the hypothetical case T — 0, when the excitations are 
completely localized, the three-particle states do not con- 
tribute to the Raman intensity. Thus the T term plays 
an important role providing both a finite bandwidth and 
coupling to higher bound state modes. Let us also men- 
tion that three-particle bound states have been predicted 
theoretically in quasi-lD systems |ll[ . Our work presents 
the first experimental evidence for their existence in a 
(2D) dimcrized spin system. 

The bound state modes experience strong damping 
with increasing temperature (see Fig. 1(b)), due to scat- 
tering on thermally excited triplet states. The decrease of 
the scattering intensity Ib with temperature is governed 

— A 

by the density of excitations: Ib(T) oc (1 — Ae k B T ) , 
where A is a constant. In Fig. 2(a) the integrated inten- 
sity of the bound states is shown together with a fit based 
on the equation above (A=215.7) and fixing A=34 K 
from experimental data. The macroscopic temperature- 
induced population of triplet states has a direct effect 
on the excitation spectrum. For intermediate temper- 
atures (5-100 K, typical spectrum shown at 7 K in 
Fig. 1(a)) a broad continuum of scattering is observed 
that is of magnetic origin ^M . However, in comparison 
with a 2D antifcrromagnet with a maximum of the two- 
magnon continuum at F, max —2.8 J jl9| its energy is very 
small. The reason for the maximum appearing at only 
50 cm -1 «2A is the smallness of the one-particle band- 
width in SrCu2 (1303)2- The temperature dependence of 
the two-magnon scattering intensity shows a pronounced 
increase in the temperature range where the intensity of 
the bound states is dropping (see the open triangles in 
Fig. 2(b)). At higher temperatures (not shown here) it 
decreases again and is then superimposed on the inten- 
sity of quasielastic scattering. In Fig. 2(b) we compare 
the integrated scattering intensity with the behavior of 
the magnetic susceptibility x(T). Calculations of the lat- 
ter quantity exist, but their agreement with experiment 
is not fully satisfactory ||l]j4| and consequently we take 
x(T) directly from experiment 0. Notice that both the 
two-magnon continuum intensity and x(T) exhibit strong 
variations on the same temperature scale, even though 
we are not aware of a rigorous sum rule relating the two 
quantities. 

Quasielastic scattering connected with fluctuations of 
the magnetic energy density is found at high tempera- 
tures (T>A) in our Raman scattering experiments. It 
has the expected Lorentzian spectral function ]2(|. To 
determine the evolution as a function of temperature 
we use the hydrodynamic form of the respective corre- 
lation function |HJ , which includes the magnetic specific 
heat c m (T) and the thermal diffusion constant Dt- In 
the high temperature approximation (hLo/ksT <C 1) the 

result is f(|: I(uj) oc if J^DtS)' 1 ' where k is thc 
scattering wave-vector. A fit to this equation can then 
be used to estimate c m (T) from the integrated intensity, 
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scaled by T 2 . In Fig. 2(c) the result of this procedure is 
found in good agreement with the measured specific heat. 

In conclusion, we have shown that the low-energy ex- 
citation spectrum of the 2D compound SrCu2(B03)2 has 
a rich and complex structure. Our main result is that 
at low temperature (T<A) the spectrum contains sev- 
eral well pronounced singlet modes which arc interpreted 
as collective two-particle and novel three-particle sing- 
let bound states of strongly localized triplets. We have 
demonstrated that since quantum fluctuations are absent 
in the ground state (due to the unique dimer arrange- 
ment), the triplet-triplet interactions can lead to large 
binding energies and scattering intensities, and conse- 
quently make the collective modes observable. At inter- 
mediate (T~A) and high (T^>A) temperatures, where 
the quantum (interaction) effects are not important, the 
light scattering is dominated by a two-magnon contin- 
uum and quasielastic contributions. We relate the spec- 
tral weight evolution in these regimes to thermodynamic 
quantities. 
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FIG. 1. Raman spectra of SrCu2(BOa)2 in two scattering 
geometries. The upper inset shows Cu-dimers of 1/4 of the 
unit cell formed by Ji (full line) and J2 (dashed line). The 
respective polarization (arrows) of incident and scattered light 
is given with respect to the crystallographic axes. The lower 
inset shows spectra for B=0 (full line) and 6 Tesla (dotted 
line) in (a b' + c) polarization. An additional A s -phonon is 
marked by an arrow. 
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FIG. 2. Mapping of Raman scattering intensities (open 
symbols) on thermodynamic quantities (lines) as a function 
of temperature. The scattering intensities and the thermody- 
namic data from Ref. mm are normalized to their maxima. 
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